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Abstract 



• A finite rank separable approximation for the quasiparticle RPA with Skyrme 



interactions is applied to study the low lying quadrupole and octupole states 
in some S isotopes and giant resonances in some spherical nuclei. It is shown 
that characteristics calculated within the suggested approach are in a good 
agreement with available experimental data. 

I. INTRODUCTION 

The random phase approximation (RPA) with the self-consistent mean-field 

derived with making use of the Gogny's interaction or the Skyrme-type interactions 
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is nowadays one of the standard tool to perform nuclear structure calculations. Many 
properties of the nuclear collective states can be described successfully within such models 

Due to the anharmonicity of vibrations there is a coupling between one-phonon and 
more complex states The main difficulty is that the complexity of calculations beyond 
the standard RPA increases rapidly with the size of the configuration space and one has 
to work within limited spaces. It is well known that using simple separable forces one can 
perform calculations of nuclear characteristics in very large configuration spaces since there 
is no need to diagonalize matrices whose dimensions grow with the size of configuration 
space. For example, the well-known quasiparticle-phonon model (QPM) belongs to such 
a model. Very detailed predictions can be made by QPM for nuclei away from closed shells 



That is why a finite rank approximation for the particle-hole (p-h) interaction resulting 
from Skyrme-type forces has been suggested in our previous work |T^. Thus, the self- 
consistent mean field can be calculated in the standard way with the original Skyrme in- 
teraction whereas the RPA solutions would be obtained with the finite rank approximation 
to the p-h matrix elements. It was found that the finite rank approximation reproduces 
reasonably well the dipole and quadrupole strength distributions in Ar isotopes fTSll . 

Recently, we extended the finite rank approximation for p-h interactions of Skyrme type 
to take into account pairing [|19|. We tested our approach to calculate characteristics of the 
low-lying quadrupole and octupole states in some spherical nuclei. In this paper we apply 
our approach to study the low lying quadrupole and octupole states in some S isotopes. 
Choosing as examples some spherical nuclei we demonstrate an ability of the method to 
describe correctly the strength distributions in a broad excitation energy interval. 
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II. BASIC FORMULAE AND DETAILS OF CALCULATIONS 



We start from the effective Skyrme interaction and use the notation of Ref. | 20[ 



containing exphcit density dependence and all spin-exchange terms. The single-particle 
spectrum is calculated within the HF method. The continuous part of the single-particle 
spectrum is discretized by diagonalizing the HF hamiltonian on the harmonic oscillator basis 
H\. The p-h residual interaction Vres corresponding to the Skyrme force and including both 



direct and exchange terms can be obtained as the second derivative of the energy density 



functional with respect to the density |2^. Following our previous papers ||T8| , p!9[| we simplify 
Vres by approximating it by its Landau-Migdal form. For Skyrme interactions all Landau 
parameters Fi,Gi, F^, G'l with / > 1 are zero. Here, we keep only the / = terms in Vres and 
in the coordinate representation one can write it in the following form: 

Vresin, ra) = ATg-i [Fo(ri) + G'o(ri)(Tia2 + (Fo(ri) + Go(ri)cTi(T2)rir2] <5(ri - rs) (1) 

where cTj and are the spin and isospin operators, and A^o = 2kpm* /n'^fi'^ with kp 
and m* standing for the Fermi momentum and nucleon effective mass. The expressions for 
Fo,Go, Fq,G'q in terms of the Skyrme force parameters can be found in Ref. [1^]. Because 
of the density dependence of the interaction the Landau parameters of Eq.(|I|) are functions 
of the coordinate r. 

The p-h residual interaction can be presented as a sum of N separable terms. To illustrate 
a procedure for making the finite rank approximation we examine only the contribution of 
the term Fq. In what follows we use the second quantized representation and Kes can be 
written as: 

Vres = ^ H ^1234 : afa^a^a^ : (2) 

^ 1234 

where af (ai) is the particle creation (annihilation) operator and 1 denotes the quantum 
numbers (ni/ijirrii), 

Vi234 = J 0i(ri)0;(r2)Kes(ri,r2)03(ri)04(r2)dridr2, (3) 



^234 = E^"'0'lll^^lb3)(j2||lj||j4)/(J1J2J3J4) X (4) 
JM 



In the above equation, (ji||yj||j3) is the reduced matrix element of the spherical har- 
monics Fj^, J = \J2J + 1, and I[^i\j23z3^ is the radial integral: 

where u(r) is the radial part of the HF single-particle wavefunction. As it is shown in 
T^,|T9[ the radial integrals can be calculated accurately by choosing a large enough cutoff 



radius R and using a A^-point integration Gauss formula with abscissas and weights Wk- 
Hjihisji) ^ Nq^^ Y1 ^^^"'^^^^ Uj, {rk)uj,, (rfc)Mj3 {rk)uj^ (r^) (6) 



^ fc=l 'k 

So we employ the hamiltonian including an average nuclear HF field, pairing interactions, 
the isoscalar and isovector particle-hole residual forces in the finite rank separable form 
||T9| . This hamiltonian has the same form as the QPM hamiltonian with separable terms 
PJ23[| , but in contrast to the QPM all parameters of this hamiltonian are expressed through 
parameters of the Skyrme forces. 

In what follows we work in the quasiparticle representation defined by the canonical 
Bogoliubov transformation: 



a 



Ujttf^ + {-ly ""Vjaj^rn- (7) 



The single-particle states are specified by the quantum numbers {jm) The quasiparticle 
energies, the chemical potentials, the energy gap and the coefficients u,v of the Bogoliubov 
transformations are determined from the BCS equations. 
We introduce the phonon creation operators 

QU = ^ E {X^f AHjj'; A/i) - {-r}'-^Y^ A{j/; X - /.)) (8) 
jj' 



where 
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(9) 



The index A denotes total angular momentum and fi is its z-projection in the laboratory 
system. One assumes that the QRPA ground state is the phonon vacuum | 0), 
i-e. Qx^i I 0) =0. We define the excited states for this approximation by Qx^i \ 0). 

Making use of the linearized equation-of-motion approach one can derive the QRPA 
equations 

( 



A B 
B -A 





(x] 




(x) 






= W 




I 









(10) 



In QRPA problems there appear two types of interaction matrix elements, the matrix 
related to forward-going graphs '^(ji/)(j2i^) ^^'^ matrix related to backward-going graphs 
^ijij')U2j'2)' Solutions of this set of hnear equations yield the eigen-energies and the ampli- 
tudes X, Y of the excited states. A dimension of the matrixes i3 is a space size of the 
two-quasiparticle configurations. Expressions for A, B and X, Y are given in [|1^] . 

Using the finite rank approximation we need to invert a matrix having a dimension 
AN X AN independently of the configuration space size. One can find a prescription how to 
solve the system (|TUp within our approach in |TB|,|TP|]. The QRPA equations in the QPM 
have the same form as the equations derived within our approach [ll8| , [T9|| , but the single- 
particle spectrum and parameters of the p-h residual interaction are calculated making use 
of the Skyrme forces. 



In this work we use the standard parametrization SIII ||2J] of the Skyrme force. Spherical 
symmetry is assumed for the HF ground states. It is well known ^^^^ constant gap 

approximation leads to an overestimating of occupation probabilities for subshells that are 
far from the Fermi level and it is necessary to introduce a cut-off in the single-particle space. 
Above this cut-off subshells don't participate in the pairing effect. In our calculations we 
choose the BCS subspace to include all subshells lying below 5 MeV. The pairing constants 
are fixed to reproduce the odd-even mass difference of neighboring nuclei. In order to perform 
RPA calculations, the single-particle continuum is discretized by diagonalizing the HF 
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hamiltonian on a basis of twelve harmonic oscillator shells and cutting off the single-particle 
spectra at the energy of 160 MeV. This is sufficient to exhaust practically all the energy- 



weighted sum rule. Our investigations |T9| enable us to conclude that A^=45 is enough for 
multipolarities A < 3 in nuclei with A < 208. Increasing A^, for example, up to A^=60 in 
208p|^ changes results for energies and transition probabilities not more than by 1%, so all 
calculations in what follows have been done with A^=45. Our calculations show that, for 
the normal parity states one can neglect the spin-multipole interactions as a rule and this 
reduces by a factor 2 the total matrix dimension. For example, for the octupole excitations 
in 206p]^ UYg| need to invert a matrix having a dimension 2N=90 instead of diagonalizing 



a 1376 X 1376 matrix as it would be without the finite rank approximation. For light nuclei 
the reduction of matrix dimensions due to the finite rank approximation is 3 or 4. So, for 
heavy nuclei our approach gives a large gain in comparison with an exact diagonalization. 
It is worth to point out that after solving the RPA problem with a separable interaction, 
to take into account the coupling with two-phonon configurations demands to diagonalize a 
matrix having a size that does not exceed 40 for the giant resonance calculations in heavy 
nuclei whereas one would need to diagonalize a matrix with a dimension of the order of a 
few thousands at least for a non-separable case. 

III. RESULTS OF CALCULATIONS 

As a first example we examine the 2i and 3j~ state energies and transition probabilities 
in some S isotopes. The results of our calculations for the energies and B(E2)-values and 
the experimental data [^] are shown in Table 1. One can see that there is a rather good 



agreement with experimental data. Results of our calculations for S isotopes are close to 
those of QRPA with Skyrme forces [Q. The evolution of the B(E2)-values in the S isotopes 
demonstrates clearly the pairing effects. The experimental and calculated B(E2)- values in 
^^S are two times less than those in ^'^'^'^S. The neutron shell closure leads to the vanishing 
of the neutron pairing and a reduction of the proton gap. As a result there is a remarkable 
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reduction of the E2 transition probability in ^^S. Some overestimate of the energies in ^^'^^S 
indicates that there is room for two-phonon effects. The study of the influence of two-phonon 
configurations on properties of the low-lying states within our approach is in progress now. 

Results of our calculations for the 3j~ energies and the transition probabilities B(E3) 
are compared with experimental data ET^ in Table 2. Generally there is a good agreement 



between theory and experiment. 

An additional information about the structure of the first 2^, 3~ states can be extracted 
by looking at the ratio of the multipole transition matrix elements M„/Mp that depend on 
the relative contributions of the proton and neutron configurations. In the framework of 
the collective model for isoscalar excitations this ratio is equal to M„/Mp = N/Z and any 
deviation from this value can indicate an isovector character of the state. The Mn/Mp ratio 
can be determined experimentally by using different external probes [^,^,0. Recently 



26| , QRPA calculations of the Mn/Mp ratios for the 2i states in some S isotopes have 
been done. The predicted results are in good agreement with experimental data [^]. Our 
calculated values of the Mn/Mp ratios for the 2f and 3j~ states are shown in Tables 1 and 2, 



respectively. Our results support the conclusions of Refs. EH] about the isovector character 



of the 2t states in ^^S. As one can see from Table 2 our calculations predict that the Mn/Mp 
ratios for the 3j~ states are rather close to N/Z, thus indicating their isoscalar character. 

To test our approach for high lying states we examine the dipole strength distributions. 
The calculated dipole strength distributions (GDR) in ^^Ar, ^^^Sn and ^°^Pb are displayed in 
Fig. 1. For the energy centroids (mi/mo) we get 19.9 MeV, 15.8 MeV and 12.7 MeV in ^^Ar, 
^^^Sn and ^''^Pb respectively. The calculated energy centroid for ^°^Pb is in a satisfactory 
agreement with the experimental value |^1| (13.4 MeV). The values of energy centroids for 
^^Ar and ^^^Sn are rather close to the empirical systematics p2[ Ec = 31. 2y4~^/'^ + 20.6^4"^/^ 



MeV. For ^^Ar the QRPA gives results that are very similar to our previous calculations with 
the particle-hole RPA ||TB[ because the influence of pairing on the giant resonance properties 



is weak. It is worth to mention that experimental data for the giant resonances in light 
nuclei are very scarce. 



The octupole strength distribution in ^°^Pb is rather well studied in many experiments 
33| , p^ . The calculated octupole strength distribution up to the excitation energy 35 MeV is 



shown in Fig. 2. According to experimental data [^] for the 3i state in ^osp]^ ^^le excitation 
energy equals to = 2.62 MeV and the energy-weighted sum rule (EWSR) is exhausted by 
20.4% that can be compared with the calculated values = 2.66 MeV and EWSR=21%. 
For the low-energy octupole resonance below 7.5 MeV our calculation gives the centroid 
energy Ec = 5.96MeV and EWSR=12% and experimental values are 5.4 MeV and 15.2% 
accordingly. For the high-energy octupole resonance we get values Ec = 20.9 MeV and 
EWSR=61% that are in a good agreement with experimental findings Ec = 20.5 ± 1 MeV 



and EWSR=75 ± 15% One can conclude that present calculations reproduce correctly 
not only the 3f characteristics, but the whole octupole strength distribution in ^°^Pb. 



IV. CONCLUSION 

A finite rank separable approximation for the QRPA calculations with Skyrme interac- 
tions that was proposed in our previous work is applied to study the evolution of dipole, 
quadrupole and octupole excitations in several nuclei. It is shown that the suggested ap- 
proach enables one to reduce remarkably the dimensions of the matrices that must be in- 
verted to perform structure calculations in very large configuration spaces. 

As an illustration of the method we have calculated the energies and transition probabil- 
ities of the 1~, 2+ and 3~ states in some S, Ar, Sn and Pb isotopes. The calculated values 
are very close to those that were calculated in QRPA with the full Skyrme interactions. 
They are in an agreement with available experimental data. 
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TABLES 

TABLE I. Energies, B(E2)-values and {Mn/Mp)/{N/Z) ratios for up-transitions to the first 2+ 
states 



Nucleus 




Energy 


B(E2t) 




{Mn/Mp)/{N/Z) 






(MeV) 


(e^fm^) 








Exp. 


Theory 


Exp. 


Theory 


Exp. Theory 


32s 


2.23 


3.34 


300±13 


340 


0.94±0.16 0.92 


34s 


2.13 


2.48 


212±12 


290 


0.85±0.23 0.87 


36s 


3.29 


2.33 


104±28 


130 


0.65±0.18 0.40 


38s 


1.29 


1.55 


235±30 


300 


1.09±0.29 0.73 


TABLE IL 


Ener 


gies, B(E3)-values and (M„/Mp)/(Ar/Z) ratios for up-transitions to the first 


3~ states 












Nucleus 




Energy 


B(E3t) 




(M„/Mp)/(iV/Z) 






(MeV) 


(e^fm^) 








Exp. 


Theory 


Exp. 


Theory 


Theory 


32s 


5.01 


7.37 


12700±2000 


8900 


0.89 


34s 


4.62 


5.66 


8000±2000 


8500 


1.06 


36s 


4.19 


3.86 


8000±3000 


7200 


1.15 


38s 




5.68 




6200 


1.01 
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FIGURES 

FIG. 1. Strength distributions of the GDR in ^^Ar, i^^Sn and ^ospb 
FIG. 2. The octupole strength distribution in 208p|-, 
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